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Abstract 



We discuss the functional representation of fermions, and obtain exact 
expressions for wave-functionals of the Schwinger model. Known features of 
the model such as bosonization and the vacuum angle arise naturally. Con- 
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l/"") | trary to expectations, the vacuum wave-functional does not simplify at large 

distances, but it may be reconstructed as a large time limit of the Schrodinger 
■ functional, which has an expansion in local terms. The functional Schrodinger 

equation reduces to a set of algebraic equations for the coefficients of these 
terms. These ideas generalize to a numerical approach to QCD in higher 
dimensions. 
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43 ■ I. INTRODUCTION 

Although the Schrodinger representation is a natural context for developing non- 
perturbative methods in field theory, and provides useful analogs of quantum mechanical 
techniques and concepts, it has not received the attention it deserves. In particular, the 
existence of the Schrodinger equation was not shown until the work of Symanzik [|IJ . How- 
ever there has been a growing interest in the subject, as field theorists attempt to relate 
increasingly sophisticated non-perturbative ideas to the more concrete calculations which 
may be performed, for example, in lattice gauge theory. The Schrodinger representation 
also arises naturally in the study of field theories on space-times with boundaries 0, and 
has recently been applied to investigation of the AdS/CFT correspondence ||. 

Symanzik's work was restricted to perturbation theory, but in principle the Schrodinger 
equation can be solved non-perturbatively. In [ffl, it was shown how to reduce the 
Schrodinger equation for scalar fields to an infinite set of coupled algebraic equations. Wave- 
functionals (WF's) are non-local in general, but in the presence of a mass gap they undergo 
a simplification; for slowly-varying fields they have a derivative expansion, the terms of 
which are local expressions involving the fields and their derivatives at a single space-time 
point. By exploiting analyticity in a complex scale parameter, it was shown in H how this 
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expansion could be used to reconstruct physics at all length scales. In this paper we extend 
these ideas to fermionic fields. 

Since the physical spectrum of the Schwinger model has a mass gap, we calculated 
its vacuum WF in the expectation that it would have a local expansion of the form just 
described. However, we will find that the vacuum WF contains non-physical massless modes 
which enforce gauge invariance by screening certain large-distance configurations; because of 
this, the vacuum WF cannot be expanded in local terms. This failure of the local expansion 
property seems to be a general feature of gauge-invariant theories; for example a similar result 
was found for Yang- Mills in |J. In spite of this, however, we can still obtain the vacuum 
WF from a local expansion by re-expressing it as the large-time limit of the Schrodinger 
WF. The latter has a local expansion for all field theories, even those containing massless 
particles, because its time parameter acts as an inverse mass for all propagating modes. We 
will verify this explicitly for the Schwinger model. 

The Schwinger model has received extensive study over the past several decades, largely 
because it illustrates many features of QCD, such as confinement and chiral symmetry break- 
ing, in an understandable way. It also provides a useful model for testing new techniques. 

This paper is organized as follows. In section 2 we discuss the Schrodinger representation 
for fermions coupled to gauge fields. Extending the discussion of the fermion representation 
given by Floreanini and Jackiw in [[?]], we argue that WF's depend on fields which are 
constrained by a local chiral- type projection implemented by operators Q±. We identify all 
possible choices of Q± in arbitrary dimension, arguing that they are parametrized by a single 
complex number. We find a path-integral expression for the vacuum WF of a gauge theory 
coupled to fermions. In section 3 we explicitly obtain the vacuum WF of the Schwinger 
model in two ways, by solving the path-integral, and by solving the Schrodinger equation. 
The path- integral solution provides a particularly straightforward derivation of this result. 
We discuss various physical issues such as the vacuum angle, bosonization, confinement, and 
vacuum condensates. In section 4 we obtain the Schrodinger functional, verify that it has 
an expansion in local terms, and show that its large time limit reproduces the result of the 
preceding section. In section 5 we present our conclusions and discuss the significance of 
these results in the context of a proposed numerical approach to higher dimensional QCD. 



Wave-functionals for fermions coupled to a gauge field may be represented as overlaps 
with a field state. We work in the Weyl gauge A = 0, so that the arguments of WF's are 
fields (Grassman- valued for fermions) defined on a surface of constant time, and WF's are 
invariant under time-independent gauge transformations of these fields. So for the vacuum 
WF, for example, we have 



where we have suppressed vector, spinor, colour and flavour indices. We will find that u and 
u are actually constrained fields, subject to a local chiral-type projection. The field state 
is chosen as follows: it is diagonal in the gauge field 



II. REPRESENTING FERMION FIELDS 



V[u,u\A] = (u,u\A\0) 



(1) 



(u,u\A\A = A(u,u\A\ 



(2) 
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and the conjugate electric field is represented by functional differentiation 

{u,u\A\E=-i—{u,u\A\. (3) 

For fermions we use the Florenanini-Jackiw representation 

1 5 
((u,u\A\4> = -j=(u + j^)((u,u\A\ 

((u,u\A\^ = -L(«t + l-)(( u ,u\A\. (4) 

The notation of indicates that u and are temporarily treated as unconstrained 
fields, eg. for the purpose of taking functional derivatives. It is obvious that @) is reducible, 
but it is important to note that its reducibility may be removed by imposing the constraints 
Q_u = v)Q + = 0, for some arbitrary projection operators Q± = |(1 ± Q), Q 2 = 1. In the 
presence of these constraints, the representation (f|) corresponds to diagonalizing Q + ip and 

(u,u\A\Q + ip = V2Q + u(u,u\A\ 

(u,u\A\ftQ- = y/2v!Q-(u,v!,A\, (5) 
and representing the other projections by functional differentiation 

(u,u\A\Q-ip = -^=Q^-^(u,u\A\ 

(u,v!,A\ftQ+ = -^=—Q+(u,u^A\. (6) 

We will sometimes use (§]) for convenience, but this is equivalent to using ([|) and 
So physical WF's will always depend only on the constrained fields, and in our final results 
the constraints will always be assumed. 

We can make the field dependence of the field state more explicit by writing[] 

(u,u\A\ =e- Tr{u ' Qu) ((u,u\A\ = (Q\exvTr[V2(u*Q^-$Q+u)+iA-E], (7) 

Note the relative factor of e _Tr( " tQu ) between (u, u\ A\ and ((u, u\ A\, which disappears 
when the constraints are imposed. (Q\ is defined by 

(Q\A=(Q\Q+$ = (Q\ftQ- = 0. (8) 

Similarly, defining \u,u',A) = exp Tr[y/2(vti> _ fo u ) + iA ■ E]\Q) with \Q) satisfying 
A\Q) = Q-4>\Q) = $Q+\Q) = 0, we have 



In this section traces are functional, ie. over spatial indices as well as colour and flavour, if 
present. 
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i>\u,u\A)) = -j={u - j^)\u,u\A)) 

ft\u,u\A)) = i(«t - l-)\ u ,u\A)). (9) 

Next we wish to define an inner-product. The definitions (0) and (||) give rise to the 
following equations: 

= {4> ] - V2u^Q+((u,u\A\v,v\A}) 
= Q + $ + V2v)((u,v*,A\v,v\A)) 
= (ft + y/2rf)Q- ((it, u\ A\v, v\ A)) 

= Q44>-V2u)((u,ut,A\v,vt,A)), (10) 

where we have used the canonical commutation relations. The field operators in (|T0"| ) may 
be represented by either pair of fields, and the equations solved to give 

((it, it f , A\v, v\ A)) = exp Tr(u f Qu - u ] 2Q_v + v ] 2Q + u + v ] Qv)5{A - A). (11) 

Hence 

(it, u\ A\v, v\ A) = exp Tr(2w f u - 2u ] v)5{A - A), (12) 

which differs slightly from the inner-product given in and elsewhere. To calculate 
inner-products of WF's we use the partition of unity 

1 = J DADuDu* 1 DvDv^e Tr ( 2vtu ~ 2utv ^ \v, v\ A) (u, u\ A\ . (13) 

Now for many purposes it is convenient to require physical WF's to be gauge invariant, 
so we wish to choose Q in such a way that this is made possible. In general, we could 
consider Q to have some gauge field dependence. However, any field dependence, or non- 
locality of Q will cause the representation to transform non-trivially under time-independent 
gauge transformations (this is easily seen from (|7|) and (§)). The only gauge invariant 
representations are when we take Q to be a local, field independent operator. 

Several authors, inspired by the similarity of Q with a choice of Dirac sea, have taken Q± 
to be (non-local) projections P± onto +ve/-ve energy eigenstates. The resulting WF's are 
invariant under time-independent gauge transformations of the fields, but in general they do 
not satisfy Gauss' law |7]|§. However, if we take Q to satisfy the conditions above, Gauss' 
law is always automatically satisfied. 

We can obtain further restrictions on the choice of Q by requiring that the overlap (1) is 
well-defined and non-vanishing; since |it, u\ A) is not a physical state this is not guaranteed 
to be the case. Consider the vacuum. With P± defined as above, it satisfies 

ftP+\0) = P-$\0) = 0. (14) 
The conditions (§) and ( |T4T ) and the representation (Q) lead to the following equations: 
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These have the solution 



P-(iu + —^)((u,u\A\0) = 
(«t + A)p + ((„, ut„A|o> = 
Q+(u-—^)((u,u\A\0) = 

(v!--^)q-((u,u\a\o) = o. (15) 



(It, It 



t^lo) = (g|fie utM "ex P Tr(U-i)|0), (16) 



where M = A~ X C = G~ X A, and we define A = {P_,Q + }, C = [P-,Q+]. The constant 
of integration B corresponds to the determinant of the Dirac operator, and solving ( |l"5"D is 
equivalent to performing the fermion integration, since (|K|) no longer involves fermion field 
operators. We assume that A^ 1 exists; we can always choose Q± so that this is the case. 

We can reproduce this result in a path integral representation. The vacuum is represented 
as lim^oo e~ Ht \S), where S is any physical state not orthogonal to the vacuum. This leads, 
via the usual correspondence of matrix elements with path integrals, to the expression 

\P [it, u\ A] = J DAiD^Di/;exp^-S B - S F 

+Tr[v / 2(M t Q_V> - ^Q + u) + u*Qu + iA ■ A] |. (17) 

Here Sb and Sp are the bosonic and fermionic parts of the Euclidean action, and the integral 
is evaluated with the following boundary conditions implied by @: 

A i \ t=0 = Q + ^\ t=0 = ^Q-\t= = 0. (18) 

The easiest way to impose the boundary conditions is by point splitting the equal time 
propagators which arise in (|TTD . We choose for ip^Q+ and Q-ip to be propagated forward in 
time, and ip^Q-, Q+ip backwards (the direction of propagation of the gauge field is arbitrary). 
This is effected by adding boundary terms 

{ip ] _ e Q + ip) + iTr(A" e • A) (19) 

to the action, where ip_ e (x,t) = ip^(x,t — e), etc. and we let e — > 0. After adding these 
terms, we are free to choose any boundary conditions we like, or to integrate freely over the 
boundary values; the result is equivalent to evaluating ( |TT| ) with the boundary conditions 
implied by @. 

The fermion integration is now easily done by standard methods; we obtain the path 
integral version of ( |16D 

^[u,u\A}= J DA t detD exp^-SB + TT^A^Cu + tiA-A)^]}. (20) 

D is the Dirac operator. We can also write down the vacuum functional for free fermions or 
fermions in a classical background field: 
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= expTrfwU^Cw 



(21) 



We will now investigate the possible choices for Q. The solution we found above depended 
on the existence of the operator A~ 1 C; for local Q if this does not exist, (JX3T) has no non- 
singular solution. We can write A~ 1 C — (Q + P)~ 1 P(Q + P), so we need Q + P to be 
invertible for all momenta. All Q satisfying this condition are given in Appendix E of [|TJ; 
they take the form 



Q = cry + z67°7 5 + C7 5 
where a 2 + b 2 + c 2 = 1, with the additional conditions 

(l-b 2 Y l (-a±ibc) £ [l,oo), 
0. 



(22) 



(23) 



for m > 0, and b 2 ^ 1 for m 
We can rewrite fl22| ) as 



1 



Q + = (a 1 + b 1 1 )-(l + a' 1 ° + b' 1 1 ), 



(24) 



Q- = \{\-a'<f -6V)(a 7 ° + &7 1 ), (25) 

where a' = and 6' = ^qf§> from (||) and (||) we see that this is equivalent to choosing 
Q± = \{1 ±a' 7 ° ife'7 1 ). Since a' 2 + b' 2 = 1, this means that we can set c = in (|22"D without 
loss of generality. Thus the possible choices for Q are parametrized by the single, arbitrary 
complex number z = a — ib (which could in principle be taken to have local variations in 
space) . 

Specializing to two dimensions we can take 7 ° = a , 7 1 = a 2 , 7 5 = cr 3 , in which case 

± = i(J_, f). (26) 

In two dimensions, the general form of a physical WF is 

°° 1 a [' ■ f Xn A 

^ = (u,u\A\^ f ) = J2~ II dx n dy n ^(x n )-? 5 u(y n )e ieJ y™ f a {x 1 ,y 1 , . . . ,x a ,y a ), (27) 

a=0 a - n=l J 

where the Wilson lines e* e ^« guarantee gauge invariance and f a are arbitrary functions. 
It is useful to note that in the presence of the constraints on u and w', two dimensional 
gamma matrices have only one effective degree of freedom; ie. in (|27|) we could have written 
7 or 7 1 in place of 7 5 . 

Using this and (|13|) we can calculate (dropping a divergent constant) 

00 1 r 

(V g \Vf) =g°f° + Y,- I ^x^ye^g^x^yi,, . . . ^^y^YTix^y^ . . . ,x a ,y a ). (28) 

a=l ®" 



6 



III. THE SCHWINGER MODEL VACUUM 



A. Path Integral Solution 



In (1+1) dimensions gauge theories are sufficiently simple to admit analytical treatment, 
and provide a great deal of information about the general structure of WF's. Consider the 
Schwinger model, with Euclidean action 



J d 2 x^ 1 ■ (d + ieA)^ + -F^F^). 



(29) 



We take gamma matrices to be Euclidean, 7 M , 7^ = 5^. For massless fermions in (1+1) 
dimensions we can obtain the fermion determinant by integrating the anomaly ||; with a 
U(l) gauge group the result is 



det D = exp 



2tt 



(30) 



where = d^t] + e^dycj). P± are obtained as equal time limits of the Dirac propagator 

P± = lim D~ l , (31) 

and 

D _ gie(»?+7 5 9i)^0^^ e -ie(??+7 5 (/)) 

At t = we want A\ = A, which implies that r\ = f x (A + <f)). For we choose the boundary 
condition 4>\t=o = 0, corresponding to the vanishing of the vacuum angle, as we will explain 
later. Thus we have 

P± = I Ux -y)± ^P^y) exp{^e [{A + 0)}. (33) 



u 



Here V denotes the principal part. The operator A 1 C is easily found; fl2HD becomes 
^[u,u\A] = J £>0expj— ^(d 4 + M 2 d 2 )(j) 



+ — I dxdy 



h 5 u(y)V—^—exp\ie[ (A + <fi)\ 
x-y 1 Jy J 



(34) 



where M 2 = — . We can integrate this to give 



oo i a 



a=0 °- n=0 



2 f 1 

- / dx n dy n u\x n )rfu{y n )V 

7i J x n - y n 



ie f X A 

e J y 



where 



x exp j J2 H x i -Vi)- J2 " x i) + ~ Vi) 

»J=1 j>i=l 



(35) 



(36) 



It may be explicitly verified that ( |31)D satisfies both Gauss' law and the Schrodinger equation. 
Note that this expression is both UV and IR convergent; apart from a subtraction of zero- 
point energy, no renormalization is necessary in this model. 
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B. The Vacuum Angle 



The WF (RH) is strictly gauge-invariant, so it corresponds to the vacuum with vanishing 



ie9 
2^ 



vacuum angle. However, a non-zero vacuum angle may be obtained by inserting a term 

dxe^F^ (37) 

into the path- integral expression (|TT|). This term is proportional to the "instanton number", 
which is not quantized, as a result of the boundary conditions. We have e^F^ = —2d 2 <f), 
but on the other hand A = — /"^ dtd 2 (f), so that 

J d 2 xe iiV F llv = 2 J dxA. (38) 

Hence 

V g [u,u\A] = y[u,u\A}e^rf dxA (39) 

so that under a large gauge transformation J dxA — ► J dx^4 + — , we have *# — > \E , ee 2me '. 
This phase is twice what we might expect; because of the relation ( |38|) large gauge transfor- 
mations change the instanton number by multiples of two. To change the instanton number 
by one unit, we must perform a large gauge transformation of winding-number one half 
/ dxA — > / dxA+ -. These "half- integer" transformations were first discussed in |10| ) but 



should not be thought of as contributing to the vacuum degeneracy JT~TJ] . The periodicity of 
9 is explained by the fact that it is energetically favourable to create an electron-positron 
pair when 9 > 2n ||12|| . 



Because of the chiral anomaly, there is another representation for the theta vacuum; by 
performing a chiral rotation 

^ e^ 5/ V, ^ -> ^ t e-^ 5/2 , (40) 

in (|17D the term ( ^7|) is cancelled. This returns us to the original expression, but with Q 
replaced by Q' = Qe t8 ^ and u, chirally rotated, ie. 

■$ e [u,u\A] = y[e- ie ^/ 2 u,v}e ie ^/ 2 ,A}. (41) 

This implies that 

(^16(^,^)1*1) = (^^(e^ 5 / 2 ^,^^-^ 5 / 2 )!* 2 ), (42) 

a fact which is useful for calculating expectation values in the ^-vacuum. 

All of these results are unaffected by giving the fermions a non-zero mass, except that 
the mass term is also altered by the chiral rotation (|40"|); m — > me ldl . This indicates that 
the vacuum angle has a true physical significance in this case, whereas for massless fermions 
it can be altered by a suitable redefinition of the fields. 
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C. Solution from the Schrodinger Equation 



A solution for the vacuum WF of the Schwinger model (on a circle) was given in [13|, 
where it was found by solving the Schrodinger equation. We now wish to show how a similar 
calculation can reproduce the result of the last section. Define 

J M (x,y) = ft(x)M$(y)e ie Jy A . (43) 

Gauss' law implies that d x E(x) ~ ejo( x ) on physical states, where jo = ^tp- By Fourier 
transforming, this in turn gives us 

J dxE 2 {x) ~ j dxdyj {x)] {y)V{x -y)+F 2 , (44) 



where F = J dxE and 



v ^ = rJp"- (45) 



Thus the Hamiltonian acting on physical states can be written as 

H — J dxdyjo(x)jo(y)V(x — y) — i J dx lim d y J y 5(x, y) + F 2 . (46) 

Note that the fermion interaction is completely determined by Gauss' law. The F term is 
represented by a ^' d dxA ^ , an d may be diagonalized by inserting a factor of e~ J dxA into the 

WF's, ie. it gives the vacuum angle. For the moment we will take this to vanish, so that 
this term may be ignored. 

Gauge invariant WF's have the form (p7[). If the distributions f a are taken to be trans- 
lation invariant, then ( p7|) is simultaneously a zero eigenstate of momentum. Inspired by 
the path integral calculation of the last section we make the Ansatz 

/ a (xi,2/i, . . .,x a ,y a ) 

a [ a a i 

= 2 J] A(x n - y n ) exp I J - y t ) - Yl [ $ (^ ~ x i) + ~ Vi)\ \ > ( 47 ) 

n=l [«,i=l J>t=l J 

where A(x) = ^V-^. The equal-time limits of the free-fermion propagator are |(1 ±7 5 A), 
and $ is the two-point function of the effective theory which results from integrating out the 
fermions . In a more general theory we would include higher n-point functions in (4"7D, but 



there are no higher n-point functions in this CcLSG, clS the effective theory is non-interacting. 
We take $ symmetric and $(0) = 0, so we can write 

$(x) = J dpC(p)(l —cos(px)), (48) 

and we have only to determine C(p). Let us write 

= J2 II J dx n dy n u j (x n )'j 5 u(y n )e eJ y^ f a (x 1 ,y 1 , . . .,x a ,y a ). (49) 



a=0 ' n=l ' 
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We will interpret sums and products from 1 to as and 1 respectively. An explicit 
computation, as in [13|, gives 

f a (xi,y h ...,x a , y a ) = V(x h y x , . . . , x a , y a )f a (x u y 1} . . . , x a , y a ) 

-i / dx lim 9 y / a+1 (xi, y x , ...,x a , y a , y, x) 

J y-*x 



/a 
dxlimd y Y<f a+1 (xi,yi 

V ^ X 6=1 



%bi y, yb, ■ ■ ■ , x a} y a ) 



-i f a 1 (?i,yi,..., ftb, jib,---, x a , y a )5'(x b - y h ), 

b=l 

where ft denotes the omission of that argument, and we define 

a a 

V(xi, yi,..., X a , Va) = J2 V ( X i ~ Vi) ~ YI Wi X i ~ X o) + V (Vi ~ Vj)\ > 

i,j=l j>i=l 

and V(x) = 2(V(0) -V(x)). 

We wish to verify that \l/ solves the Schrodinger equation 



(50) 



(51) 



For a = (|52|) reduces to 

E = -i 



- x)e* {y ~ x) 



)} 



dx l™{ d y( A (y 

and gives us the divergent vacuum energy. For a = 1 we find 

£ A(xi - yi) = V(xx - yi)A(xx - y x ) 

+i J dxA(xi — x)A{x — yi)(&(x — y\) — $'(x — x\)) 

-iA{x x -Vx) J dx lim d y (A(y - x ) e s(wi,2/i)^ ; 



(52) 



(53) 



(54) 



where g(x,y,xi,yi) = $(?/ — x) + $(xi — x) + $(y — yi) — $(y — x-y) — $(x — yi), and we 
have used / dxA(xi — x)A(x — x-z) = 8{x\ — a^)- Now by Taylor expanding a test function 
f(x) it is clear that 



limA(i)/(i) = -/'(0) + ... 



71 

— I 



limA'(,)/(x) = -/"(0) + 

In 



(55) 



The ellipses represent potentially divergent terms which contribute to the vacuum energy 
but cancel in (|54|). Hence we have 



— 1 



J dxlimdy (A{y - x)e 9{x ' y > xl ' Vl) ) = E + — J dx{&{x - y x ) - &{x - x^) 2 , (56) 



and 
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— J dx(&(x - vx) - &{x - x 1 )) 2 = J dpp 2 C{p) 2 {l - cos(p(x 1 - Vl ))). (57) 
Now the identity 

dxA(x 1 - x)A(x - yi)f(x) = 6(x! - 

+A( Xl - yi ) J ^ c (p)/(p)( e -*« _ e -w) (58) 

corresponds in momentum space to the identity e(p)e(p + g) = 1 + e(p)e(g) — e(q)e(p + q) 
(e is the step function). This implies that 



2 / dxA(x\ — x)A(x — yi)($'(x — yi) — — xi)) 

= 2A( Xl - y x ) J dp\p\C{p){l - cos(p(x 1 - Vl ))). (59) 
Thus (0) reduces to the quadratic equation 



e 



2 



Tip 1 



+ 2\p\C(p)-p'C( P y = 0, (60) 



and to get a normalizable WF we must take the appropriate root 

m = lpl ' ^ TW . (en 

This reproduces the result (55). It is straightforward to show that fl52|) is satisfied by this 
solution for all a > 1; using the identities (l57f) and (|58|) the a = n equation may be reduced 
to multiple copies of the a = 1 equation. 



D. Bosonization 

The confining nature of the theory is easily seen from (|35"D; since the potential (|45|) goes 
like |x| at large distances, configurations in which a field source goes off to infinity without 
an accompanying anti-source, are exponentially damped. Infrared slavery and asymptotic 
freedom are similarly seen by scaling the momentum in (fi5D; the effect of such scaling is to 
make the coupling constant large for small momenta and vice-versa. 

Bosonization may be understood by transfering the sources to currents instead of fields. 
By a gauge transformation, we can always set ^4 = 0. Now consider the following object: 

[f] = y Du j DueI d2 ' xu ^ 5u ^[u, 0]. (62) 

Inserting the path-integral expression (^4|) and performing the u, integrations gives (with 
boundary conditions as before) 

[ D<j>Difi ■ J D^ e -s s -^+v^(^ t 7 5 V0 . ( 63 ) 
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This is a functional whose argument couples to a local fermion current. Using ( |30|) the 
fermion integration yields 



where (p = e<p — ^(j) an d 4> — ^\^d\J^ with J = 0, Ji = i^/n^ft^). Hence <fi = i^^{^5{t)) 
and performing the remaining integration, we find 

y[£] = e -*f d3 *V = ®W*t, (65) 

which is the vacuum WF of a free boson field of mass M. 

From this we can easily reproduce the standard correspondence between bosonic and 
fermionic operators. ( j63| ) implies that 



1 5 



>,/,t r>M,U^i„-SB-s F +^(^t^) 



= J D<j)D^Dilj^j 5 ^e 

= J Du ( Duef d2xuf ^ u ^j 5 ^[u, 0]. (66) 

The fermion operators are represented as in ([5]) and @. Now define 

[f , x] = / D^I^itye -5 * -5 ^^^^"^^. (67) 

Performing the fermion integration we find J\ = iy/Tr£5(t), as before, and Jo = y/nx5(t), so 
that = ^|(^(t)) + 0f|( x 5(t)). Thus 

= e Jdx(-^V-d^+M^-x^~h 9 i(-9 2 +M 2 r 1/2 x)_ (68) 

Hence 



| T^lW^^'V^K ^"^ t , 0]. (69) 



If we represent £ as — then fl6"6|) and flBTf ) imply that 



DutlW^^YV^fa, — Cw f , 0] 



/ WlW^^^L^d^K 0]. (70) 

j \/7r 2 



It follows that ip'y fj "ijj ~ ^Le^d^ on all physical states. This equivalence may be exploited 
in a number of ways; for example it immediately allows us to identify the creation and 
annihilation operators of the theory 



12 



a±(p) = i (^ 1/2 J (P)/P T ^ 1/2 j 5 (p)) , (71) 

where jo(p) and jsQo) are the Fourier transforms of jo(x) = ffi(x)ip(x) and js(x) = 
■0"l"(:r)7 5 ?/>(:r) respectively. Also, in calculating physical expectation values 

(^10(^)1^) (72) 

the operator O can be represented in terms of equivalent bosonic operators; the resulting 
calculations are generally much easier. To illustrate this, consider the chiral condensate 

W V)e (6\6) 

(73) 



(0|0) ' 

where we have used (pE2"|). Now the bosonic operators corresponding to the chiral densities 

ip{l ± 7 5 )V' are 

^f- : e ±{ ^ : (74) 

(this can be verified in a similar way to the other bosonization formulae). The normal- 
ordering is perfomed with respect to the scale M, and 7 is the Euler constant. Thus x/je 101 -0 
corresponds to Mjf- : cos(v^4~7r£ + 6) : whose vacuum expectation value is easily ascertained 
from (|65|). We have 

^^) e = ^-cos9, (75) 



which is the well-known result 14 



IV. THE LOCAL EXPANSION 
A. The Schrodinger Functional 

We will now describe how the result ( pop can be reconstructed from a derivative expan- 
sion, whose terms are local expressions in the fields. This allows the techniques of fU|, 
which are generalizable to higher dimensions, to be used to provide an alternative solution 
of the theory. Naively, the presence of a mass gap in the theory should mean that all prop- 
agators are exponentially damped at large distances, so that the logarithm of (p5|) should 
reduce for slowly varying fields to a sum of local terms, ie. integrals over finite powers of 
the fields and their derivatives, evaluated at a single spatial point. Unfortunately, because 
of gauge-invariance (^) contains massless modes not appearing in the physical spectrum, 
and this simplification of the vacuum WF does not occur. This may be seen by noting the 
existence of screened large-distance configurations. 

We will get around this problem by considering instead the Schrodinger functional 
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V T [u,ifi,u,vt] = (u,u ] \e~ Hr \u,u ] ). (76) 

The vacuum functional is just the r — > oo limit of this object. We can also extract excited 
states by inserting a basis of energy eigenstates 

tf T [0, ~Y,*E[u,v!]e- ET . (77) 

The parameter r, which corresponds to Euclidean time, acts as an inverse mass for all states 
of the theory — physical or otherwise. Provided we work with fields that vary slowly on the 
scale of r, the logarithm of ty T has an expansion in positive powers of r, each term of which 
is itself a finite sum of local terms. As we will see, the large r behaviour can be reconstructed 
from this expansion. 

We will begin by finding the free fermion solution. The Schrodinger equation is 

- JU T = HV T , (78) 

where 

H=]-Jd d x^ + j^)h(u + ^), h = i^d x . (79) 

We also have the initial condition 

^o = (Q\u,v?)) = eI dx{u ' Qu \ (80) 



If we put 



then (|50D becomes 



^ = e / < fa(«tr(T)«) j ( gl ) 



f = -I(i-r>(i + r), (82) 

with T(O) = Q. This equation is solved by 

r = (£ + Q_)(£-Q_)-\ (83) 



and substitution into (52) gives £ = hT,, S(0) = Q + , with the solution 

S = e-^ Tdx Q + . (84) 

Hence we obtain 

r = Q + 2Q_7 5 tanh(-ir^)g+, (85) 

which has a derivative expansion for small r, as promised. It is also easy to find the (non- 
local) large-time limit: 
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• 5 1 

lim T(x-y) = Q + 2Q^V Q+, (86) 

r ^°° 7T x — y 

which coincides with the solution T = A _1 C that we found before. 

Now consider the interacting theory. Substituting the result fl50"D into the time-dependent 
Schrodinger equation flTRD gives 

- f a (xi,yi, ■ ■ ■ , x a , y a ) = -Eof a (xi,yi, ...,x a , y a ) + V(x x ,y x , ...,x a , y a )f a (xi,y 1 , ...,x a , y a ) 

-i / dxlimd y f a+1 (x 1 ,y 1 , . . . ,x a ,y a ,y,x) 

/a 
dx lim d y V f a+1 (xi,yx,..., x b , x, y,y b ,..., x a , y a ) 

6=1 

a 

-*Z)/ a_1 ( x i'^i' • • -ifaitib, ■ ■ .,x a ,y a )5'(x b - y b ), (87) 

6=1 

Here we have made H regular by subtracting the zero-point energy (p3[), ie. H — > H — Eq; 
this may be achieved by normal-ordering. In order that we recover the solution fl8"5|) in the 
free-fermion limit e — > 0, we use the Ansatz (H7I) with A(x — y) = tanh(ir9 x )5(x — y). The 
a = part of (p7|) is of course trivial; for a = 1 it becomes 

= A(xi - ?/i)<9 T $(xi - yi) + V(x x - y 1 )A(x 1 - y x ) 

+i J dxA(xi — x)A(x — yi)(&(x — yx) — $'(x — Xx)) 

-iA(xx - yx) J dxlimdy (A(y - x )e?( x > v < xuVl ' i ) _ # A(a;i - Vl ). (88) 
In obtaining this we used the identity 

A(xx — yx) = i(l ~ tanh 2 (frd))5'(:ri — ?/i) 

= « 5' (xi — yx) — i J dxA(xx — x)A'(x — yx). (89) 

Now for small x 

ITT I 7T 2 X \ I 

A(x) = — cosech + 0(x), (90) 

r \ r J irx 

so that 

\imA(x)f(x) = -f(0) + ... 

limA'(x)/(x) = ^/"(0) + ... (91) 
The ellipses have the same meaning as in (0). Hence as before 

- i J dxlimdy (A(y - x)e 9{x ' y ' Xl ' yi) ) = E - J dpp 2 C(p,T) 2 (l - cos(p(xx - yx)))- (92) 
On the other hand, we find that 
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i J dxA(x x - x)A(x - yx)(<&'(x - y x ) - ®{x - x x )) = 

2A(x x - y x ) J dpcoth(rp)pC(p, r)(l - cos(p(x x -y x ))). (93) 

Thus fl88f) reduces to 

M 2 

C(p, t) + 2p coth(pr)C(p, r) - p 2 C 2 (p, r) + — = 0. (94) 
The initial condition flSOf ) is satisfied if C(p, 0) = 0. For small r we can expand C as 

oo 

C = X c„r-, (95) 

n=l 

and substituting into ( |94"|) leads to a recursion relation which is easily solved for the coeffi- 
cients c n , which are polynomials in positive powers of p, divided by l/p 2 . We note that fl9"I|) 
reduces to ( |60|) for large r, so that we recover the vacuum functional, as expected. 

Just as for the vacuum functional, it is easily shown that the solution given by ( |9^| ) 
satisfies (^) for all a. 



B. Solution from a local Ansatz 

The Ansatz ([T7]) is very useful in two dimensions, but it is not clear that it generalizes in 
any way to higher dimensions. Thus we will now describe an Ansatz which does generalize. 

Provided we work with fields that vary slowly on the scale of r, A(x) = td,nh(id x T)5(x) 
can be expanded in derivatives of the delta function. It follows from ([|7]) that the logarithm 
of \l/ r has a local expansion (this should also be true in higher dimensions). This allows us 
to write log\l/ T in the form (0), where we take the f a to be local: 

oo 

f a (p X ,...,P2a)= £ & ni P? 1 • • • b n2aP Ta 
ni,...,ri2a=0 

f a (x 1 ,y 1 ,...,X a ,y a ) = [ ^ . . . C ^5(p 1 + . . . + p 2a )f a (PU ■ ■ ■ ,P2a)- (96) 

J 2,71 271 

Since we can gauge away the electromagnetic field, this leads to the following Ansatz 

„ oo 

^ r = 1 + / dx X a^ ) r% t (x)7 5 9>(x) 

n,q=0 

„ oo 

+ dxdy X 4trV(x) 7 5 (9>(x)) M t( y ) 7 5(^ M ( 2/ )) + ..., (97) 

n,m, <j=0 

where the coefficients ajj... are obtained from the b n in the obvious way. This Ansatz can be 
inserted directly into the Schrodinger equation (|3TD, which thus reduces to an infinite set of 
algebraic equations (a finite set for each order in r). We have explicitly verified for the first 
few terms the that the local expansion which results from solving the Schrodinger equation 
in this way is the same as that which is obtained by expanding the small-r solution found 
in the last section. 
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Now the local expansion depended on expanding \I/ T in positive powers of r for small 
t, but we can reconstruct the large r behaviour from a knowledge of the local expansion 
alone. If we evaluate Hf T for scaled fields u(x/ \fp) and u\x/J~p) then it can be analytically 
continued to the complex p-plane with the negative real axis removed. The proof of this is a 
straightforward generalization of the arguments of [fH|. Then by using Cauchy's theorem, 
we can relate the value of the functional at p = 1 to the value at large p, where the fields 
are slowly varying, so that the r-expansion converges for large r. Specifically, we study the 
following integral: 

^X^T^^/^^^/Vp)]- (98) 

If we take C to be a circle of sufficiently large radius we can evaluate this integral by 
inserting the local expansion of ty T for any value of r. Alternatively, we can collapse the 
contour around the point p — 1, and the negative real axis, the latter contribution being 
exponentially suppressed for large A. Thus we have successfully expressed \l/ r for large r in 
terms of the local expansion. In practise we will wish to truncate the expansion at some 
order, in which case the truncation error is minimized by taking A to be large but finite [§J . 

The Ansatz (|96| ) generalizes straightforwardly to higher dimensions, and could form the 
starting point for a numerical approach to (3+1) dimensional QCD. 

V. DISCUSSION 

In this paper we have shown how the degeneracy of the Floreanini-Jackiw representation 
for fermions is related to a choice of diagonalization of the fermionic fields. We have used 
this to derive formal expressions for gauge-invariant wave-functionals in the Schrodinger 
representation, and obtained the vacuum wave-functional of the Schwinger model, both 
from a path-integral expression and as a solution of the Schrodinger equation. All known 
features of the model arise naturally in this context. 

In general, the wave-functionals of a theory with a mass gap should be expandable in 
local quantities at large distances. We have seen how the existence of screened large-distance 
configurations spoils this feature in the Schwinger model. However, we have obtained the 
Schrodinger wave-functional and shown that it does have an expansion in local quantities. 
This may be obtained directly from the Schrodinger equation, which reduces to a set of 
algebraic equations. An analogous technique can be applied to gauge theories in higher 
dimensions, where the Schrodinger functional should also have such an expansion. Physical 
wave-functionals may be reconstructed by considering its behaviour for large values of the 
time parameter, which acts as an infrared regulator. 

The techniques of this paper generalize in a straightforward way to other two dimensional 
theories, for example involving massive fermions or a non-abelian gauge group. Various 
perturbative and non-perturbative approximation schemes present themselves. And though 
the detailed structure of the solutions is specific to two dimensions, the solution via a local 
expansion is a technique which should work in a similar way for higher dimensional gauge 
theories. 

The main steps involved in implementing a numerical approach to higher dimensional 
QCD based on solving the Schrodinger equation for the local expansion are as follows. 
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First, to obtain the generalization of (|44j ) to non-Abelian fields; this is straightforward, 
and allows the Hamiltonian to be written in terms of gauge-invariant currents. Second, to 
write down the local ansatz which generalizes (0). The local gauge- invariant terms which 
make up the local expansion are easily identified. Renormalization may be dealt with by 
adding counterterms familiar from perturbation theory. Or, we can consider super symmetric 
theories in which they are absent. Finally, after finding the action of the Hamiltonian on the 
terms of the local expansion, the solution of the Schrodinger equation proceeds as described 
in section 4 above. 
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